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Abstract 



We discuss forward jet production data recently published by the HI 
and ZEUS collaborations at HERA. We review how several Monte-Carlo 
models compare to the data. QCD calculations based on the BFKL for- 
malism and on fixed NLO perturbation theory with and without resolved 
virtual photons are described. 
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Abstract. We discuss forward jet production data recently published by the HI 
and ZEUS collaborations at HERA. We review how several Monte-Carlo models 
compare to the data. QCD calculations based on the BFKL formalism and on fixed 
NLO perturbation theory with and without resolved virtual photons are described. 

1 Introduction 

It is by now firmly established that the proton structure function Fn(x, Q 2 ) 
shows a steep rise in the small x-Bjorken region, i.e., below x — I0~ 3 
This rise is compatible with DGLAP evolution ||, where F2 is fitted at a 
fixed input scale Qq and then evolved up to Q 2 by summing up InQ 2 terms. 
In the small x region an alternative way to describe the evolution might be 
to sum up ln(I/a;) terms, as it is done in the BFKL approach The BFKL 
equation actually directly predicts the behaviour of F2 as a function of x. 
An attempt to unify these two approaches is the CCFM evolution equation 
||, which reproduces both, the DGLAP and the BFKL behaviour, in their 
respective regimes of validity. 

The existing data on Fi do not allow to unambigously determine whether 
the BFKL mechanism is needed in the ir-range covered by HERA. Therefore, 
alternatively the cross section for forward jet production in deep inelastic 
scattering (DIS) has been proposed as a particularly sensitive means to in- 
vestigate the parton dynamics at small x Q. The proposal is based on the 
observation that the DGLAP and BFKL equations predict different ordering 
of the transverse momenta fcj^ along the parton cascade developing in DIS 
jet production, see Fig. 1. While the DGLAP equation predicts a strong kx,i 
ordering and only a weak ordering in the longitudinal momentum fractions 
Xij i.e., 

Q 2 = k\ n » . . . 3> k T1 and x = x n < ... <X\ (I) 

the BFKL approach predicts a strong ordering in the Xi, but no ordering 
in kx,i- The idea is now to observe jets at very small x, with 1 « 11 in 
the forward region with Q 2 ~ k 2 ^ 1 . In this way the DGLAP evolution is 
suppressed, whereas the BFKL evolution is left active. 

The forward jet cross section has been measured recently at HERA |?],|| 
and in the following I will review the various attempts that have been made 
to describe this data. 
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Fig. 1. Evolution of the transverse momenta kr,i and the longitudinal momentum 
fractions Xi along the ladder contributing to DIS jet production. 



2 Data and Monte-Carlo models 



The HI and ZEUS collaborations have measured forward jet cross sections 
at small x for rather similar kinematical conditions The jet selection 

criteria and kinematical cuts are summarized in Tab. 1. 

In Fig. H and || the collected data from ZEUS and HI in the forward 
region is shown together with predictions from various Monte-Carlo models. 
The data from both groups is on the hadron level. For the ZEUS data, the 
Monte-Carlo models ARIADNE §, LEPTO fo), HERWIG (n| and LDC g| 
have been used for predictions. ARIADNE includes one of the main features 
present in the BFKL approach, which is the absence of the strong hp ordering. 
LDC is based on the CCFM approach and finally LEPTO and HERWIG are 



Table 1. Forward jet selection criteria by HI and ZEUS 



HI cuts ZEUS cuts 

E' e > 11 GeV E' e > 10 GeV 

Ve > 0.1 y e > 0.1 

E T ,jet > 3.5 (5) GeV E T , ]e t > 5 GeV 

1.7 < t] jet < 2.8 t] jet < 2.6 

0.5 < E 2 T ^JQ 2 < 2 0.5 < E 2 T ^JQ 2 < 2 

x jet > 0.035 x jet > 0.036 



Forward jet production at HERA 3 



ZEUS 1995 



ZEUS 1995 




ZEUS Data 
ARIADNE 4.08 
LEPTO 6.5 (SCI) 
HERWIG 5.9 
LDC 1.0 



10' 



10 




ZEUS Data |^ 
. ARIADNE 4.08 
, LEPTO 6.5 (SCI) 

HERWIG 5.9 

LDC 1.0 



Tr- 



io 



10 



Fig. 2. Forward jet data of ZEUS M as a function of Bjorken-x for Er jet > 5 GeV. 
(a) linear scale, (b) logarithmic scale. The shaded band gives the error due to the 
uncertainty of the jet energy scale. The data is compared to various Monte-Carlo 
models. 
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Fig. 3. Forward jet data of HI || as a function of Bjorken-a; for two Er jet cuts of 
3.5 GeV and 5.0 GeV. (a) and (c) contain Monte-Carlo model predictions, whereas 
(b) and (d) show the results of a LO BFKL (full) and a fixed NLO (dashed) calcu- 
lation. 
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based on conventional leading log DGLAP evolution. Except for HERWIG, 
the same Monte-Carlo models are shown in the HI plots. Instead of HERWIG, 
the RAPGAP model |L3| is shown, which is also based on DGLAP evolution 
but contains an additional resolved virtual photon component. 

As is clear from Fig. ||, ARIADNE describes the forward jet cross section 
reasonably well, apart from the smallest x-bin where it gives slightly too small 
cross sections. The other three used models predict cross sections which lie 
significantly below the data. Similar results can be extracted from Fig. |^ 
(a) and (c). LDC and LEPTO lie below the data by a factor of 2, whereas 
ARIADNE gives a reasonable good description. As an interesting result, also 
RAPGAP gives a good description of the HI data. 



3 BFKL approach 



Of course, attempts have been made to calculate the forward jet cross section 
directly within the BFKL formalism. BFKL calculations in LO by Kwiecihsky 
et al., Bartels et al. and TangjITj] based on the BFKL approach overshoot 
the older forward jet data BM. This can also be seen in Fig. | (b) and 
(d), where the LO BFKL calculation on the parton level |l6| (full line) is 
compared to the recent HI data ||. These older calculations suffer, however, 
from several deficiencies. They are asymptotic and do not contain the correct 
kinematic constraints of the produced jets. Furthermore they do not allow the 
implementation of a jet algorithm as used in the experimental analysis. Also 
NLO ln(l/x) terms in the BFKL kernel jlTj predict large negative corrections 
which are expected to reduce the forward cross section as well. 

Recently the BFKL calculations have been improved by taking into ac- 
count higher order consistency conditions as a way of including sub-leading 
corrections to the BFKL equation The consistency constraint (CC) re- 
quires that the virtuality of the emitted gluons along the chain should arise 
predominantly from the transverse components of momentum. By including 
this CC the authors in Jji| claim to recover a dominant part of higher order 
effects. Furthermore, the CC is said to subsume energy-momentum conser- 
vation over a wide range of the allowed phase space, which is another source 
of sub-leading contributions. 

Including the CC conditions, good agreement between the predictions and 
the forward jet data is found, for both the HI and the ZEUS data, as shown 
in Fig. ||. The predictions depend on the choice of scale and on an additional 
infrared cut-off parameter fcrj. However, the kg dependence is much less than 
the uncertainty due to the choice of scales. Further details about the BFKL 
calculations from can be found in these proceedings Jl9| , also describing 
the calculation for forward ir° production. 
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Fig. 4. Forward jet data of HI fl and ZEUS |?J compared to prediction based 
on the BFKL formalism including sub-leading corrections j^]. The three curves 
correspond to three choices of scales and infrared cut-off. 



4 Fixed NLO QCD calculations 

Calculations in fixed order perturbation theory have already been performed 
for the older HI forward jet measurement ]l5| by Mirkes and Zeppenfeld p0| 
using their fixed order program MEP JET pjj . The calculations where done 
in next-to-leading order (NLO) accuracy i.e., taking the matrix elements up 
to 0(a 2 s ) into account. It was found that the NLO calculations are a factor 
of 2 to 4 below the data. Similar predictions have been made with help of the 
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DISENT program |22| for the more recent HI measurements ||, shown as the 
dashed lines in Fig. |3| (b) and (d), which confirm the earlier findings. Since 
the forward jet data can be succesfully described by the RAPGAP model 
which includes resolved virtual photons, it is fair to ask whether also fixed 
order calculations including a resolved virtual photon component will be able 
to describe the data. 



4.1 Low Q 2 jet production in NLO 

NLO calculations for jet production with slightly off-shell direct and resolved 



virtual photons have become available recently 1 23 , extending calculations 
done in the photoproduction regime [^4| . The NLO calculations |^3| are per- 
formed with the phase space slicing method. As is well known the higher 
order (in a s ) contributions to the direct and resolved cross sections have 
infrared and collinear singularities. For the real corrections singular and non- 
singular regions of phase space are separated by a technical cut-off parameter 
y s . Both, real and virtual corrections, are regularized by going to d dimen- 
sions. The NLO corrections to the direct process become singular in the limit 
Q 2 — ► in the initial state on the real photon side. For Q 2 = these pho- 
ton initial state singularities are usually also evaluated with the dimensional 
regularization method. Then the singular contributions appear as poles in 
e = (4 — d)/2 multiplied with the splitting function P ql and have the form 
— -Pqry multiplied with the LO matrix elements for quark-parton scattering. 
These singular contributions are absorbed into PDF's / a / 7 (x) of the real 
photon. For Q 2 ^ the corresponding contributions are replaced by 

- l -P ql ^-Hs/Q 2 )P ql (2) 

where ^fs is the cm. energy of the photon-parton subprocess. These terms 
are finite as long as Q 2 ^ and can be evaluated with d = 4 dimensions, but 
become large for small Q 2 , which suggests to absorb them as terms propor- 
tional to \n(M 2 / Q 2 ) in the PDF of the virtual photon. Parametrizations of 
the virtual photon have been provided by several groups ]2j| . By this absorp- 
tion the PDF of the virtual photon becomes dependent on M 7 , which is the 
factorization scale of the virtual photon, in analogy to the real photon case. 
Of course, this absorption of large terms is necessary only for Q 2 <C M 2 . In all 
other cases the direct cross section can be calculated without the subtraction 
and the additional resolved contribution. M 2 will be of the order of E\. But 
also when Q 2 ~ M 2 , one can perform this subtraction. Then the subtracted 
term will be added again in the resolved contribution, so that the sum of 
the two cross sections remains unchanged. In this way also the dependence 
of the cross section on M 2 must cancel, as long as the resolved contribution 
is calculated in LO only. 

In the general formula for the deep-inelastic scattering cross section, one 
has two contributions, the transverse (da^ b ) and the longitudinal part (da^ b ). 
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Since only the transverse part has the initial-state collinear singularity the 



subtraction in |23 has been performed only in the matrix element which 



contributes to da^ b . Therefore the longitudinal PDF's f£, are not needed. 
It is also well known that da^ b vanishes for Q 2 — > 0. The calculation of 
the resolved cross section including NLO corrections proceeds as for real 
photoproduction at Q 2 = 0, except that the cross section is calculated also 
for final state variables in the virtual photon-proton center-of-mass system. 

The NLO calculations in the low Q 2 region are implemented in the fixed 
order program JETVIP |26| . Various measurements at HERA in which the 
jet analysis has been done with JETVIP point to the presence of a resolved 
virtual photon component up to moderate virtualities of Q 2 ~ 5 GeV 2 p?)]. 



4.2 Comparison to forward jet data 

Recently, we have performed a NLO calculation including the virtual resolved 
photon for the forward jet region [^8| with the help of JETVIP. The results 
for the ZEUS kinematical conditions are shown in Fig. |5| a,b. In Fig. |5| a we 
plotted the full 0(a 2 ) inclusive two-jet cross section (DIS) as a function of 
x for three different scales = [i\ = MI 2 + Q 2 ,M 2 + Q 2 and M 2 /3 + Q 2 
with a fixed M 2 = 50 GeV 2 related to the mean E T of the forward jet 
and compared them with the measured points from ZEUS §. The choice 
Hp > Q 2 is mandatory if we want to include a resolved contribution. Similar 
to the results obtained with MEPJET and DISENT, the NLO direct cross 
section is by a factor 2 to 4 too small compared to the data. The variation 
inside the assumed range of scales is small, so that also with a reasonable 
change of scales we can not get agreement with the data. In Fig. |s| b we 
show the corresponding forward jet cross sections with the NLO resolved 
contribution included, labeled DIR5+RES, again for the three different scales 
fi as in Fig. [| a. Now we find good agreement with the ZEUS data. The scale 
dependence is not so large that we must fear our results not to be trustworthy. 

In Fig. |H| c,d we show the results compared to the HI data || obtained 
with Et > 3.5 GeV in the HERA system. In the plot on the left the data are 
compared with the pure NLO direct prediction, which turns out to be too 
small by a similar factor as observed in the comparison with the ZEUS data. 
In Fig. |^ d the forward jet cross section is plotted with the NLO resolved 
contribution included in the way described above. We find good agreement 
with the HI data inside the scale variation window M 2 /3 + Q 2 < /.i 2 < 
3M 2 + Q 2 . We have also compared the predictions with the data from the 
larger Et cut, namely Et > 5.0 GeV, and found similar good agreement p8| . 

As described in [^8) the NLO resolved contribution supplies higher order 
terms in two ways, first through the NLO corrections in the hard scattering 
cross section and second in the leading logarithmic approximation by evolv- 
ing the PDF's of the virtual photon to the chosen factorization scale. This 
way the logarithms in E T /Q 2 are summed, which, however, in the consid- 
ered kinematical region is not an important effect numerically. Therefore, the 
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Fig. 5. Dijet cross section in the forward region compared to HERA data: (a) and 
(b) ZEUS; (c) and (d) HI. (a) NLO DIS, E T > 5 GeV; (b) NLO DIR S +RES, 
E T > 5 GeV; (c) NLO DIS, E T > 3.5 GeV; (d) NLO DIR S +RES, E T > 3.5 GeV. 



enhancement of the NLO direct cross section through inclusion of resolved 
processes in NLO is mainly due to the convolution of the point-like term 
in the photon PDF with the NLO resolved matrix elements, which gives an 
approximation to the NNLO direct cross section without resolved contribu- 
tions. One can therefore speculate that the forward jet cross section could be 
described within a fixed NNLO calculation, using only direct photons. 



5 Conclusions 

We conclude that two alternative models exist for describing the forward jet 
data. The BFKL calculation including CC describes the data well. This is 
supportet by the ARIADNE model, where strong kr ordering is absent and 
which also describes the data. It is however not clear, why the LDC model 
does not decribe the data, although BFKL dynamics should be included in 
the respective regime of validity. 
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Similarly, the NLO theory with a resolved virtual photon contribution 
as an approximation of the NNLO DIS cross section, which is presently not 
available, gives a good description of the forward jet data. Assumably, the 
RAPGAP model produces these higher order effects through parton shower 
contributions in the resolved cross section. 

Acknowledgments. I thank the organizers for the kind invitation and the 
pleasant workshop athmospherc. The results of section 4 where obtained in 
collaboration with G. Kramer. 
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